arXiv:1509.00265vl [nucl-th] 1 Sep 2015 


EPJ Web of Conferences will be set by the publisher 
DOI: will be set by the publisher 

© Owned by the authors, published by EDP Sciences, 2015 


Electric dipole polarizability: from few- to many-body systems 

Mirko Miorelli^'^’®, Sonia Bacca^’^''^, Nir Barnea"^, Gaute Hagen®’®, Giuseppina Orlandini^’®, and 
Thomas Papenbrock®’® 

^ TRIUMF, 4004 Wesbrook Mall, Vancouver, British Columbia, Canada VST 2A3 

^Department of Physics and Astronomy, University of British Columbia, Vancouver, British Columbia, 
Canada VST 1Z4 

^Department of Physics and Astronomy, University of Manitoba, Winnipeg, Manitoba, Canada R3T 2N2 
‘'Racah Institute of Physics, Hebrew University, 91904 Jerusalem, Israel 
^Physics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 
^ Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 3799S, USA 
^ Dipartimento di Fisica, Universita di Trento, Via Sommarive 14, 1-38123 Trento, Italy 
^Istituto Nazionale di Fisica Nucleare, TIFPA, Via Sommarive 14, 1-38123 Trento, Italy 


Abstract. We review the Lorentz integral transform coupled-cluster method for the 
calculation of the electric dipole polarizability. We benchmark our results with exact 
hyperspherical harmonics calculations for '*He and then we move to a heavier nucleus 
studying *®0. We observe that the implemented chiral nucleon-nucleon interaction at 
next-to-next-to-next-to-leading order underestimates the electric dipole polarizability. 


1 Introduction 

The electric dipole polarizability is a fundamental quantity to understand nuclear dynamics and is 
related to the response of the nucleus to an external electric field. Typically, it is measured via photo¬ 
absorption reactions ID, Compton scattering H or hadronic processes, such as ip,p') reactions a. 
Ab-initio calculations of this observable were traditionally available only for very light nuclei. It is 
our goal to extend such studies to the medium mass regime using coupled-cluster theory. 

2 Review of the method 

The electric dipole polarizability can be calculated from the dipole response function R{a)) as an 
inverse energy weighted sum rule 



( 1 ) 


In the coupled-cluster formalism HI the dipole response function is 


R{to) = ^0'* e^|n«)<nL|e '^&e'^\0R}6{E„ - Eq - tu). 


( 2 ) 
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where 0 is the dipole excitation operator, f is the cluster operator defined in coupled-cluster theory 
a, and (0^,1, |0«) ({«/,!,!««)) are the left and right reference ground states (excited states), respectively. 
Using Eq. Q we can then rewrite the electric dipole polarizability as 

ao = 2a<Oz.|0^|'P«), (3) 

where 0 = is the similarity transformed operator and '!'« is the solution of the Schrodinger- 

like equation 

(H - AEoWk) = 0|O«), (4) 

where AEq is the correlation energy of the ground state. This equation is obtained using the Lorentz 
integral transform method © in the coupled-cluster formulation 0121. Similarly to what was done in 
HI El, the polarizability can be expressed as a continued fraction of the Lanczos coefficients cni 


ao = 2a<Oi|0^0|OR> 


ao - 


1 



(5) 


If one is able to calculate the response function of Eq. then the polarizability could be obtained 
using the definition in Eq. Q. However, using Eq. Q allows us to calculate the polarizability directly 
from the Lanczos coefficients, thus avoiding further uncertainties coming from the inversion of the 
transform. Eurthermore, due to the nature of the continued fraction, convergence is reached with less 
than 100 Lanczos steps. Eq. Q is an exact result, free of any approximation, if the cluster expansion 
is performed up to A-particle-A-hole excitations, A being the total number of nucleons. However, in 
practical calculations the cluster expansion is truncated since a full expansion is not feasible due to the 
very demanding computational cost. In the results we present here we truncated the cluster operators 
at the singles-and-doubles excitations level and dub this approximation scheme LIT-CCSD. 


3 Results and Conclusions 

3.1 Benchmark in ^He 

While the LIT-CCSD method has been already benchmarked with the effective interaction hyperspher- 
ical harmonic (EIHH) method im for the dipole response function lH, here we show a benchmark for 
the electric dipole polarizability calculated using Eq. Q with a chiral nucleon-nucleon force derived 
at next-to-next-to-next-to-next-to-leading order IIT2l . 

In Eigure we show the electric dipole polarizability of "^He as a function of the model space size 
Nmax- The green band represents the combined experimental values from insiiii which have been 
extracted by fitting to experimental photo-absorption data MM- The continuous line (brown) is the 
polarizability obtained with the EIHH method. The lines with scatter-symbols are the results from the 
LIT-CCSD method obtained for different values of the frequency hQ. of the harmonic oscillator basis 
used to expand the potential. As expected, the polarizability is independent of the oscillator frequency 
as one increases N^ax, and converges rapidly to a value ao = 0.0815 fm^, which compares fairly well 
with the EIHH © result of ao - 0.0831 fm^ and ao - 0.0822(5) fm^ obtained using the no-core- 
shell-model with the same interaction El. Moreover, the convergence with respect to the number of 
Lanczos coefficients used in Eq. (|^ is much faster than the one observed for the response function in 
141 . Already with 15 Lanczos coefficients the error with respect to the converged value is 0.3%, which 
shows the advantage of using directly Eq. 0- The small difference between the LIT-CCSD and EIHH 
results can thus be attributed to the truncation to singles-and-doubles excitations only. 
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Figure 1. The electric dipole polarizability in ''He 
as a function of the model space size for dif¬ 
ferent values of the harmonic oscillator frequency 
hn (black, red, blue). The converged value at 
Nmax = 18 is compared with the result obtained 
from the EIHH method (brown) and the experi¬ 
mental data from GMU (green band). 


Figure 2. Comparison of the electric dipole polar¬ 
izability in '"O as a function of the model space 
size N„cix for three different values of HQ = 20,24 
and 26 MeV (black, red, blue) with the experimen¬ 
tal data from Ahrens et al. (green). 


3.2 Application to 

The very nice agreement in the benchmark for the "^He polarizability suggests that we can extend the 
calculation of ao via this method to heavier nuclei. We then consider '®0, whose response function 
has been already studied previously IIHITl. In Figure]^ we study the convergence of the polarizability 
as a function of the model space size. Again, the convergence is very fast and the final value ao - 
0.461 fm^ is independent of the underlying harmonic oscillator frequency. However, the converged 
value is rather low if compared to the experimental one of = 0.585(9) m. This trend is observed 
also in heavier nuclei such as ^°Ca and '**Ca alii, and the discrepancy is even larger when the mass 
number is increased. As shown in Figure [T] the effect of the cluster expansion truncation is tiny and, 
because of the size extensive nature of coupled-cluster theory, we expect the effect to be very small in 
heavier systems as well. Discrepancies between the calculated values and the experimental data are 
mainly attributed to deficiencies of the employed Hamiltonian. 


3.3 Conclusions 

We used the LIT-CCSD method to extend the calculation of the electric dipole polarizability from 
few- to many-body nuclei using nucleon-nucleon chiral interactions. We have shown the reliability 
of the method with benchmarks in ^He where both EIHH and LIT-CCSD calculations and the experi¬ 
mental data are in agreement. Moving to heavier systems, e.g., '®0, we observe disagreement between 
calculated and measured values. This points to the inadequacy of chiral interactions at the two-body 
level in accounting for the dipole polarizability in heavier nuclei, which we also observed for the 
radii mmol. The use of newly developed Hamiltonians with three-nucleon forces CSl calibrated on 
radii of finite nuclei can help solving this problem. Work in this direction is underway. 
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